Defect and solute properties in dilute Fe-Cr-Ni austenitic alloys: an ab initio study 
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We present results of an extensive set of first principles density functional theory calculations of 
point defect formation, binding and clustering energies in austenitic Fe with dilute concentrations 
of Cr and Ni solutes. A large set of possible coUinear magnetic structures were investigated as 
appropriate reference states for austenite. We found that the antiferromagnetic single and double 
layer structures with tetragonal relaxation of the unit cell were the most suitable reference states and 
highlighted the inherent instabilities in the ferromagnetic states. Test calculations for the presence 
and influence of non-colhnear magnetism were performed but proved mostly negative. We calculate 
the vacancy formation energy to be between 1.8 and 1.95 eV. Vacancy cluster binding was initially 
weak at 0.1 eV for divacancies but rapidly increased with additional vacancies. Clusters of up to 
six vacancies were studied and a highly stable octahedral cluster and stacking fault tetrahedron 
were found with total binding energies of 2.5 and 2.3 eV respectively. The (lOO)-dumbbell was 
found to be the most stable self-interstitial with a formation energy of between 3.2 and 3.6 eV 
and was found to form strongly bound clusters, consistent with other fee metals. Pair interaction 
models were found to be capable of capturing the trends in the defect cluster binding energy data. 
Solute-solute interactions were found to be weak in general, with a maximal positive binding of 0.1 
eV found for Ni-Ni pairs and maximum repulsion found for Cr-Cr pairs of -0.1 eV. Solute cluster 
binding was found to be consistent with a pair interaction model with Ni rich clusters being the 
most stable. Solute-defect interactions were consistent with Ni and Cr being modestly oversized and 
undersized solutes respectively, which is exactly opposite to the experimentally derived size factors 
for Ni and Cr solutes in type 316 stainless steel and in the pure materials. Ni was found to bind to 
the vacancy and to the (100) dumbbell in the tensile site by 0.1 eV and was repelled from mixed 
and compressive sites. In contrast, Cr showed a preferential binding to interstitials. Barrier heights 
for vacancy mediated migration were estimated and found to be lowest for Cr, followed by Fe and 
finally Ni. Comparison of our results with those for bcc Fe showed strong similarity for pure Fe and 
no correlation with dilute Ni and Cr. 

PACS numbers: 61.82.Bg, 71.15.Mb, 75.50.Bb 



I. INTRODUCTION 

Austenitic, face-centred cubic (fee), 7-Fe based steels 
are key materials in many applications. However, mod- 
elling the basic form, 7-Fe, is challenging because it is 
metastable under the zero-temperature conditions typi- 
cally used in quantum mechanical calculation. The high 
temperature stabilisation of fee over ferromagnetic (fm) 
body centred cubic (bcc) a-Fe at 1185K is due primarily 
to the onset of paramagnetism which contributes a con- 
siderable amount of entropy in the fee form. Phonon and 
electronic entropy also contribute, but the primacy of the 
magnetic effect is underlined by the return to paramag- 
netic bcc (S-Fe at still higher temperature (1667K). 

There has been considerable experimental effort to sta- 
bilise 7-Fe at lower temperatures, primarily by epitaxial 
growth of thin films on Cu substrates (e.g. Meyerheim^) 
and by the formation of 7-Fe precipitates by heat treat- 
ment of dilute alloys of Fe in fee Cu (e.g. Tsunoda^ and 
Hincs"^ ) . For a review of earlier work see Marsman and 
Hafner— . The principal motivation behind this effort has 
been to study the structural and magnetic properties of 



7-Fe at temperatures low enough for stable magnetic or- 
dering. The study of point defects in these systems has 
not been attempted. 

First-principles calculations have proved to be a 
very reliable method of obtaining information about 
radiation-induced defects, which previously had been un- 
reliably calculated using empirical potentials. In non- 
magnetic elements such as molybdenum and vanadium it 
was shown that the pseudopotcntial plane wave method 
reproduced experimentally-inferred self-interstitial mi- 
gration energy barriers to within 0.1 eV— , giving con- 
fidence that experimentally inaccessible quantities such 
as interstitial formation energies would also be reliable. 
This is assisted by the discovery that the strain fields 
associated with interstitials are less extensive than had 
been predicted by interatomic potentialsS'^^, such that 
calculations with supercells as small as 100 atoms can 
give near-converged solutions. 

Application of density functional theory (DFT) to 
steels is of particular interest for radiation damage ap- 
plications in which high-energy defects such as self- 
interstitials are formed. Modelling commercial steels is 



a more complicated task on account of their multicom- 
ponent nature, however, there has been much progress 
in ferritic Fe (e.g. the results of the European FP6 
PERFECT project^ and references therein) and Fe-Cr 
alloysiS.^-- which showed a number of unexpected out- 
comes. In particular, isolated Cr atoms have a small 
negative heat of solution in a-Fei2,iilj in apparent conflict 
with the phase diagram— which shows a miscibility gap. 
This conundrum was resolvedi^ii^ when it was shown 
that Cr atoms in bcc Fe exhibit strong nearest neigh- 
bour repulsion, resulting from magnetic frustration, with 
weaker repulsion still present up to 6th nearest neigh- 
bour separationii. Consequently, a dilute solution of Cr 
in a-Fe has a positive heat of formation only up to at 
most 10 at. % (at OK), a result that is consistent with 
experimental results for the heat of solution and short 
range order parameter, as discussed by Bonny et al.-^, 
but was not included in the extrapolations present in 
the phase diagram. The nonlinear variation of cohesive 
energy with concentration means that determining un- 
ambiguous energies for quantities such as the binding 
of a self-interstitial to a Cr solute proved impossible in 
concentrated alloysi^ since the calculated energy had a 
complex dependence, not only on the defect, but also on 
concentration and the atomic arrangement. 

First principles studies of 7-Fe with collinear— "— and 
non-coUineari^i^ magnetism have found many distinct 
magnetically ordered and spin spiral (meta-)stable states 
lying (approximately) between 0.08 and 0.15 eV/atom 
above the bcc ferromagnetic ground state, a-Fe. How- 
ever, to date there have not been any first principles 
calculations of the paramagnetic state of 7-Fe. There 
remains some debate about whether the paramagnetism 
is best represented as itinerant or involves localised mo- 
ments on the ions, however, throughout this paper we 
will interpret our results through the localised-moment 
picture. From a purely numerical point of view there are 
also difficulties. The Kohn-Sham functional applied in 
non-magnetic density functional theory has a single min- 
imum with respect to the wavefunctions, but collinear- 
magnetic density functional theory may have up to 2^ 
minima for an N-atom supercell, corresponding to possi- 
ble permutations of the spin. In practice, most of these 
will be unstable but one can never be sure that the low- 
est energy structure has been reached. The concept of 
metastability is also slippery, since the numerical algo- 
rithms used to find the minimum electronic energy do 
not correspond to physical pathways which the material 
can follow: the very definition of metastability is then, 
to some extent, algorithm dependent. 

The Born-Oppenheimer approximation is used to de- 
couple electronic and atomic degrees of freedom. The 
status of magnetic degrees of freedom in this approxi- 
mation is debatable. One viewpoint is that, since mag- 
netism is due to electrons, the Born-Oppenheimer surface 
is the one corresponding to the magnetic state with the 
lowest energy globally. An alternative view, which we 
adopt here, is that there are many Born-Oppenheimer 



surfaces, each corresponding to a given magnetic order- 
ing. Many of the magnetic states in fee Fe are sufficiently 
metastable to make this a useful distinction. It should, 
however, be borne in mind that DFT calculations are al- 
most exclusively performed using a numerical algorithm 
which minimises the energy of the system with respect 
to the free parameters of a set of basis functions used 
to represent the wavefunction. This algorithm does not 
correspond to any physical trajectory which the electrons 
could follow (c.f. Time dependent DFT). Consequently 
"local minimum" means a minimum from which the al- 
gorithm, in our case Block Davidson^, cannot escape in 
this basis set. It does not guarantee that the spin-state 
will be a local minimum in reality. 

To date (and to the best of our knowledge) only a 
few first principles studies of solutes and impurities have 
been performed in austenitic Fo^-'^i^^'^^' '— and only the 
work of Nazaro\s22 includes defect calculations (in order 
to study vacancy- hydrogen interactions). There exists no 
comprehensive study of point defects and their interac- 
tions in pure fee Fe or in dilute solid solutions of Ni and 
Cr in fee Fe. We present here just such a study using a 
set of magnetically ordered states to represent the para- 
magnetic state of 7-Fe. By considering more than one 
magnetic state we are able to estimate the accuracy of 
this assumption. We also intend this work to provide a 
basis for understanding the complex interactions present 
in concentrated Fe-Cr-Ni fee alloys. 

The layout of the paper is as follows. Section |TT] con- 
tains the computational details of the DFT calculations 
performed in this work. In Section IIIII we present and 
discuss our results on the bulk properties and stability 
of the magnetic reference states. The energetics of point 
defects in pure Fe and their interactions are discussed 
in Section IIVI as is their tendency to form small defect 
clusters relevant in the nucleation of microscopic defects 
such as voids and dislocation loops. The results of di- 
lute Ni and Cr solute calculations in defect free Fe and 
in interaction with point defects are presented and dis- 
cussed in Section |Vl Finally we make our conclusions in 
Section rvn 



II. COMPUTATIONAL DETAILS 

All of the following calculations have been performed 
using the mainstream DFT code VASB^ii^, a plane wave 
code that implements the projector augmented wave 
(PAW) method^? 'i?.4^ Standard PAW potentials supplied 
with VASP were used with exchange and correlation in 
the generalised gradient approximation described by the 
parametrisation of Perdew and Wang24 and spin inter- 
polation of the correlation potential provided by the im- 
proved Vosko-Wilk-Nusair scheme^. Potentials with 8, 
6 and 10 valence electrons were used for Fe, Cr and Ni 
respectively. 

The local magnetic moments on atoms were initialised 
to impose the magnetic state ordering and were then al- 



lowed to relax. The relaxed local magnetic moments were 
determined by integrating the spin density within spheres 
centred on the atoms. Sphere radii of 1.302, 1.323 and 
1.286 A were used for Fe, Cr and Ni respectively. 

Calculations of the bulk properties for Fe, Cr and Ni 
were performed with a sufficiently high plane wave cutoff 
energy (400 eV) and sufficiently dense k-point sampling 
of the Brillouin zone (e.g. 16^ Monkhorst-Pack grid for a 
conventional bcc cell) to ensure convergence of the energy 
of the system to less than ImeV per atom. 

Elastic constants were determined numerically using 
stress tensor measurements after applying small (1-2%) 
strain deformations to the equilibrium structures. These 
were cross-checked, where possible, with comparable de- 
terminations using the energy. Ensuring convergence of 
the elastic constants placed significantly higher demands 
on our calculations than for other bulk properties. A 
plane wave cutoff energy of 600 eV was used and sufficient 
k-point sampling to ensure that the elastic constants were 
converged to ± 2 GPa, e.g. a 20'^ Monkhorst-Pack grid for 
a conventional unit cell of the fee fm-LS magnetic state. 
It was also necessary to reduce the energy convergence 
criteria for electronic minimisation to 10^^ eV. 

Formation energies of defects and solutes were calcu- 
lated in supercells of 256 (±1, ±2,...) atoms, with su- 
percell dimensions held fixed at their equilibrium values 
and ionic positions free to relax. For these supercells a 
2^ k-point Monkhorst-Pack grid was used to sample the 
Brillouin zone. Test calculations showed this sampling to 
be sufficient to converge formation energies to less than 
0.03 eV in all calculations except those involving inter- 
stitial defects, where the uncertainty could be as high as 
0.07 eV. Formation energy differences and binding ener- 
gies were found to be converged to less than 0.01 eV ex- 
cept in calculations involving multiple interstitial defects 
where the error was found to be 0.03 eV. These errors are 
sufficient for our purposes. Performing calculations in a 
fixed supercell of volume, V , results in a residual pres- 
sure, P, for which an Eshelby-type elastic correction for 
finite size to the system energy^i^ of —P^V/2B, where B 
is the bulk modulus of the bulk material, can be applied 
and also serves to indicate the likely error. The size of 
these corrections is negligible in much of the work pre- 
sented here, being smaller than other sources of error in 
our calculations. However, in the largest of our vacancy 
cluster calculation and in those containing overcoordi- 
nated defects these finite size corrections are of relevance 
and are discussed in the corresponding text. 

Our data is the result of the merging of two sets of 
data, one calculated in the single layer anti-ferromagnetic 
(afml) face centred tetragonal (fct) state, the other 
mainly in the double layer anti-ferromagnetic (afmD) 
state in either of fee and fct structures but also including 
some fct ferromagnetic high spin (fm-HS) calculations. 
The plane wave cutoff energy was taken to be 300 eV for 
the afml calculations and 350 eV in all others. The non- 
convergence error in the formation energies from either of 
these plane wave cutoff energies was found to be smaller 



than the k-point sampling error. First order {N = 1) 
Methfessel and Paxton smearing-2^ of the Fermi surface 
was used throughout with smearing width, a, set to 0.3 
eV for the afml calculations and 0.2 eV in all others. 
Structural relaxations were considered converged when 
the forces on all atoms were less than 0.03 eV/A for afml 
calculations and less than 0.01 eV/A for other magnetic 
states. Formation energy calculations in the afml state 
have been performed at the equilibrium lattice parame- 
ters determined with those settings i.e. a = 3.4252A and 
c = 3. 648 A and not with those presented in TABLE H] 
The resulting differences in formation energy are, how- 
ever, negligible. 

Throughout this paper we define the formation energy, 
^'fd^x}), of a configuration containing nx atoms for 
each element, X, relative to a set of reference states for 
each element using 



Ef{{nx}) ^ Ei{nx}) -J2nxE^^ 



(1) 



where E{{nx}) is the calculated energy of the configura- 
tion and E'x ^ i^ the reference state energy for element, X. 
Here we take the reference energies to be the energies per 
atom for the pure materials i.e. Ni in its fee fm ground 
state, Cr in its bcc anti-ferromagnetic (afm) state and Fe 
in the specific ground state for the magnetic ordering we 
are studying. 

We define the binding energy between a set of n 
species, {Ai}, where a species can be a defect, solute, 
cluster of defects and solutes etc., as 



Eb{Ai,...,Ar, 



'E{{A,)-Ei{Ai,...,An) (2) 



i=l 



where E{{Ai) is the formation energy of a configuration 
containing the single species, Ai, and Ef{Ai, ...,An) is 
the formation energy of a configuration containing all 
of the species. An energetically favoured configuration 
therefore has a positive binding energy. 



III. BULK PROPERTIES AND REFERENCE 

STATES 

Density functional theory has a number of energy min- 
ima corresponding to different magnetic states with crys- 
tal structures close enough to fee to be plausible as ref- 
erence states for austenite. We take fee to mean that the 
unit cell has a = 6 = c, in contrast with fct where the 
unit cell relaxed tetragonally i.e. c ^ a — b. Our cal- 
culated energy versus volume curves for many of these 
coUinear magnetic structures are shown in FIG. [T] 

The results are consistent with, and extend, previous 
work^ii^"— 1^ and show that there are many competing 
magnetic structures very close in energy. Indeed, the en- 
ergy difference between the non-magnetic state and the 
most stable magnetic state is only 0.062 eV/atom. Test 
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FIG. 1: Energies of possible reference states from 
collinear magnetic calculations for austenitic Fe (open sym- 
bols) with ferromagnetic (fm, black, circles), non-magnetic 
(nm, red, squares) and antiferromagnetic (001) single layer 
(afml, green, upward triangles), double layer (afmD, blue, 
downward triangles) and triple layer (afmT, purple, right fac- 
ing triangles) orderings. Also shown is a magnetic (001) 
layered structure with a spin up, spin up, spin down order- 
ing (uud, beige, left facing triangles) and a magnetic ordering 
formed by taking a 4-atom fee unit cell with three spin up 
atoms and one spin down atom (fm-flip,cyan, diamonds). We 
distinguish between fee (solid curves) and fet (dashed curves) 
for the same magnetic structure. Comparison is made with 
bcc ferromagnetic Fe (black solid circles). 



calculations in both the fee afml and fee afmD struc- 
tures showed that a single flipped magnetic moment was 
(meta)-stable and cost 0.03 and 0.05 eV respectively, 
again indicating how close different magnetic structures 
are in energy. Moment flips in the fee fm-HS state proved 
costly at 0.5 eV and destabilised the fee fm-LS state with 
partial relaxation towards one of the afm states. No 
(meta-)stable moment flips were found for any of the fet 
structures. 

It is worth drawing the reader's attention to the fact 
that the afml, afmD and fm-HS magnetic states in cu- 
bic cells have been found to be unstable with respect 
to tetragonal distortion. Both afml and afmD transform 
spontaneously when the constraint is removed but the fee 
fm-HS state is an unstable equilibrium position. Tetrag- 
onal distortion away from the perfect fee structure ei- 
ther resulted in full relaxation to the ferromagnetic bcc 
ground state via the Bain path2£, i.e. by setting c/a < 1, 
or to the fet fm-HS state with c/a > 1. For atomic vol- 
umes below about 11. 4 A this unstable equilibrium posi- 
tion continuously becomes a stable equilibrium position 
for the fee fm-LS state. Overall the fet parameter space 
exhibits many local ferromagnetic minima, the most sta- 
ble being the bcc fm ground state (for c/a — l/-\/2), a 
result most succinctly presented in the contour plots of 
Spisak and Hafner-iS.. 

The elastic moduli, {-By}, of the lowest energy struc- 



tures representative of the distinct magnetic states con- 
sidered here are presented along with the lattice param- 
eters, magnetic moments and energy differences with re- 
spect to the lowest energy fet afmD state in TABLE H] 
These elastic moduli correspond to the derivatives of the 
stress tensor with respect to strains and are only equal to 
the conventional elastic moduli, {Cij}, which are propor- 
tional to the second derivative of the energy with respect 
to strain, in the case of a state of zero stress^. So for 
the fee afmD state, which is not in equilibrium, the {Bij} 
differ from the {Cy}. For all others they are identical. 
In the case of zero stress the stability of the state with 
respect to tensile and shear strains can be examined by 
calculating the eigenvalues of the elastic constant ma- 
trix. Only if these eigenvalues are positive is the state 
an energy minimum with respect to strain. For an fet 
structure the eigenvalues, {Xi}, are as follows, 

Ai = Cii — Ci2 

2A2 = (Cii + C12 + C33) 

+ ((^11+^12-^733)^ + 8^23)'/' 

2A3 = (Cii + C12 + C33) 



((Cii + C12 — C[ 



33 J 



8C 



,1/2 



13; 



as well as A4 = C44, which is doubly degenerate and 
Ae = Cee- The first three are given in TABLE H] 



Structure fee afmD fet afmD fet afml fet fm-HS 


a/A 


3.527 


3.447 


3.423 


3.418 


c/k 


3.527 


3.750 


3.658 


4.017 


c/a 


1.000 


1.088 


1.069 


1.175 


AE/eY 


0.020 


0.000 


0.014 


0.031 


fl/flB 


1.80 


1.99 


1.50 


2.40 


Bii/GPa 


224 


212 


333 


131 


Bi2/GPa 


147 


211 


241 


267 


Bi3/GPa 


81 


92 


103 


106 


Bsi/GPa 


90 


92 


103 


106 


Bss/GPa 


119 


210 


250 


289 


B44/GPa 


87 


73 


173 


56 


Bee/GPa 


108 


203 


251 


165 


Ai/GPa 


- 


1 


92 


-136 


Aa/GPa 


- 


485 


630 


503 


Aa/GPa 


- 


148 


194 


184 



TABLE I: Lattice parameters, a and c, energy per atom rel- 
ative to the fet afmD state, A_B, magnetic moment, /x, elastic 
moduli in Voight notation, Bij, and the first three eigenvalues 
of the matrix of elastic moduli, A^, for four distinct magnet- 
ically ordered structures at their energy minima. We esti- 
mate the uncertainties in the elastic moduli to be of the order 
of a few percent. The non-zero stresses present in the fee 
afmD state are, in Voight notation, cti = (72 = 3 GPa and 
as = —6 GPa. We measure the energy per atom for the fet 
afmD state to be 0.077 eV higher than the bcc fm ground 
state. 

The stability criterion for Ai, i.e. Cn > C12, is not sat- 



isfied by the fct fm-HS state, showing that it is unstable 
with respect to an orthorhombic distortion breaking the 
a = b symmetry in the lattice parameters. Despite being 
unstable it is still an equilibrium structure i.e. a saddle 
point in the energy landscape, which after a small or- 
thorhombic symmetry breaking perturbation relaxes di- 
rectly to the bcc fm ground state. Note that although 
this barrier-free double shear route from fee to bcc (via 
fct) is not the generally-considered Bain path, the rela- 
tionship between the initial and final states is identical. 
One final point worth mentioning here is that the lattice 
parameter, c — 4.017A for the fct fm-HS state is almost 
exactly v2 times the bcc fm lattice parameter, which is 
2.831 A with the settings used here. In the Bain transfor- 
mation from fee to bcc two of the lattice parameters must 
increase to exactly this value. The conclusion is that if 
we only allow one lattice parameter to increase, i.e. by 
constraining to an fct cell with c/a > 1, the system will 
still relax that lattice parameter all the way towards bcc 
and not to an intermediate value. 

The fct afmD state is stable, but is very soft with re- 
spect to a further orthorhombic distortion i.e. upon ap- 
plying the eigen-strain associated with Ai. By apply- 
ing such a strain we confirmed the existence of a wide 
minimum about zero strain with direct measurements of 
Ai yielding values between 5 and 10 GPa. This direct 
determination is still small but larger than the indirect 
measurement. 

As well as considerations of structural stability any 
prospective reference state suitable for our purposes 
should be stable with respect to the introduction of sim- 
ple point defects and solutes. The fee fm-LS state, being 
unstable even to the introduction of a vacancy, is ruled 
out, as is the fee afml state which was found to disorder 
in some calculations. The fee fm-HS state was unsur- 
prisingly found to be unstable to defects breaking the 
fee crystal symmetry. Using the fct fm-HS state, as will 
be shown, improved the situation, although it was still 
unstable to defects generating orthorhombic distortions. 
All other structures were found to be stable with respect 
to the introduction of defects. 

The final consideration in our choice of reference states 
is the usefulness of our results as a solid foundation for 
understanding the complex results in concentrated Fe- 
Cr-Ni austenitic alloys. The afml magnetic ordering has 
been found to be the only stable state for a concentrated 
austenitic alloy with composition Fe70Cr20Nil041. In 
addition there is evidence for the stability of a ferromag- 
netic state at other concentrations'*^, indicating that the 
ferromagnetic state may be stabilised by alloying. 

Given the considerations discussed above, we have con- 
centrated on four main reference states for austenitic Fe 
and dilute Cr and Ni alloys in this work. These are as 
follows: 

• The fct afml state. It exhibits stability both in 
structure and against the introduction of defects 
and is of direct relevance in the study of concen- 
trated alloys. 



• The fct afmD state. It shows reasonable structural 
stability, stability against the introduction of point 
defects and is the structure with the lowest energy 
of all collinear magnetic structures found in this 
work and elsewhere. 

• The fct fm-HS state. Despite its instabilities it 
best represents the ferromagnetic state and insights 
gleaned from its study have relevance for concen- 
trated alloy systems. It also has the closest volume 
per atom to the paramagnetic state. 

• The fee afmD state: It is stable against the intro- 
duction of point defects and is the lowest energy fee 
state. Comparison with the fct afniD state allows 
the effect of tetragonal distortion on formation and 
binding energies to be studied. 

It is worth reiterating that none of these structures 
can, by themselves, represent high temperature austenite 
or concentrated austenitic alloys. All the results in this 
paper must be accepted as approximations and only when 
a feature is common against multiple reference states are 
we confident that it is generalisable. 

Before discussing our point defect and solute calcula- 
tions in these reference states we comment on the effect 
of non-coUinear magnetism, which is known to be present 
in fee Feii2^. We calculated seven systems while allow- 
ing non-eoUinear solutions. The systems were a vacancy, 
(001) dumbbell and an octahedral interstitial in pure Fe, 
single Cr and Ni solutes, a mixed Cr-Ni (001) dumb- 
bell interstitial and a Ni solute next to a vacancy. The 
initial moments were set in an afml configuration, with 
most of the initial moments aligned coUinearly. The di- 
rections of moments on a sufficient number of atoms were 
changed, however, to perturb the system away from the 
collinear solution. All moments were free to rotate into 
non-coUinear directions in the calculations. 

Six out of these seven calculations converged to 
collinear solutions. Only in the case of a Ni solute sit- 
uated next to a vacancy were clearly non-parallel direc- 
tions of moments on atoms observed. Even then, the 
energy of the non-coUinear solution was only marginally 
lower than for the collinear solution and by an amount 
below other sources of error. Therefore, while non- 
collinearity was observed, it can be omitted for our pur- 
poses. 



IV. POINT DEFECTS IN PURE FE 

The basic quantities upon which all microstructural 
radiation-induced damage and segregation effects depend 
are the energies associated with point defect formation, 
interaction and diffusion i.e. the behaviour of the primary 
damage in radiation-induced displacement cascades. We 
calculate point defect formation and binding energies in 
pure 7-Fe based on our four different magnetic refer- 
ence states for austenite. It should be borne in mind 



that the particular magnetic ordering and the presence 
of tetragonal distortion in these states lowers their sym- 
metry relative to the perfect fee crystal structure. As 
a result, some defects which would have been uniquely 
defined in fee have multiple configurations with differ- 
ent energies. We refer to these differences as symmetry 
breaking effects and they should primarily be taken as 
an additional source of uncertainty in our calculations 
when conclusions about austenite and austenitic alloys 
are made. In order to aid the discussion of our results 
we will refer to the planes of constant moment in the 
bulk afm structures and the planes perpendicular to the 
direction of tetragonal distortion in the fct fm-HS state 
by the term, magnetic planes. The results of our point 
defect calculations are shown in TABLE ITTl 




Defect 


fee afmD 


fct afmD fct afml fct fm-HS 


Vacancy 


1.672 


1.819 


1.953 


1.692 


Octa (1) 


rlx (7) 


rlx (7) 


4.353 


3.620 


Tetra uu (2) 


3.581 


3.864 


N/A 


3.039 


Tetra ud (3) 


3.332 


3.663 


4.322 


N/A 


[110] crowdion (4) 


rlx (3) 


rlx (3) 


4.799 


3.305 


[Oil] crowdion uu (5) 


3.771 


4.255 


N/A 


(-1.791) 


[Oil] crowdion ud (6) 


3.874 


4.168 


4.818 


N/A 


[100] dumbbell (7) 


2.978 


3.316 


3.531 


(-3.181) 


[001] dumbbell (8) 


2.790 


3.195 


3.615 


2.416 


[001] dumbbell (8b) 


N/A 


N/A 


N/A 


(-2.529) 


[110] dumbbell (9) 


4.290 


4.322 


4.803 


3.288 


[Oil] dumbbell (10) 


rlx (7) 


rlx (7) 


rlx (7) 


rlx (8b) 


[111] dumbbell (11) 


rlx (3) 


rlx (3) 


4.559 


1.919 



TABLE II: Point defect formation energies, E{, in eV for the 
different (magnetic) structures studied in this work. Inter- 



stitial defect sites (numbers 1 to 6) are shown in FIG. 2(a) 



Dumbbells are identified by their axis direction in the coor- 
dinate system shown in the same figure e.g. a [001] dumbbell 
has its axis along z and perpendicular to the magnetic planes. 
Where the defect was found to be unstable the defect formed 
in the relaxation is identified. The negative formation en- 
ergies for fct fm-HS (shown in brackets) are unphysical and 
stem from the instability of this reference state with respect 
to orthorhombic distortions. The positive formation energy 
for the [001] dumbbell (for fct fm-HS) corresponds to a sym- 
metrical starting position. The negative formation energy 
(identified as 8b) resulted from the relaxation of the [Oil] 
dumbbell, which initially broke the x — y symmetry then re- 
laxed to [001]. Eshelby corrections to the formation energies 
are less than 0.07eV in all cases and at most 0.01 eV for for- 
mation energy differences. 



Problems with the fct fm-HS state are immediately ap- 
parent, i.e. negative interstitial formation energies. These 
are associated with extensive reconstruction throughout 
the unit cell. All defects exhibit significantly lower forma- 
tion energies than the other reference states. Formation 
energies in the fee afmD state are systematically lower 
than those in fct afmD by around 10%. The total con- 
figuration energies in fee are, however, still significantly 



FIG. 2: (a) (Self-)Interstitial defect positions (black circles) 
in the fcc/fct afmD magnetically ordered structure, (b) A-B 
interactions between on-site defects and solutes in the fcc/fct 
afmD magnetically ordered structure. Species A is shown in 
black, species B in grey and the surrounding Fe atoms in 
white. Arrows indicate local moments in both figures and the 
magnetic planes are shown to aid visualisation. The afmD 
state has been shown here to uniquely identify all of its dis- 
tinct defect configurations. The higher symmetry afml and 
fm-HS states share this set of configurations although many 
(e.g. tetra uu and tetra ud) will be symmetry equivalent. 



higher than the fct state. We associate the formation 
energy reductions in the fct fm-HS and fee afmD states 
with the fact that these states are not minima with re- 
spect to strain. By contrast, the fct afml and afmD states 
are minima and the influence of the defects is confined 
to the first few neighbour shells. The associated energies 
can therefore be regarded as attributable to the defect. 
We thus base our predictions for paramagnetic austenite 
on these two reference states, with their difference giving 
some indication of the error. 



A. Vacancy formation 

We find that the vacancy formation energy lies in a 
range between 1.82 and 1.95 eV. This is slightly higher 
than the typical value of one third of the cohesive energy 
observed in other transition metals^"—. The local influ- 
ence of the vacancy on the lattice was found to be highly 
dependent on the reference state. Displacements of up 
to 0.24 A towards the vacancy were found for atoms in 
the first nearest neighbour (Inn) shell in the fct fm-HS 
state whereas no displacements exceeded 0.02 A in the fct 
afml state. An intermediate value of 0.09 A was found in 
the fct afmD state. Magnetic moments were typically en- 
hanced in the Inn shell, the effect being strongest within 
a magnetic plane. Increases of 0.2 hb were found for 
the fct afmD and fm-HS states and up to 0.35 /zb in the 
afml state. At 2nn, moments within a magnetic plane 
were consistently found to be reduced in magnitude but 



by less than the enhancement seen in the Inn shell. The 
overall enhancement of moments suggests that there may 
be a significant magnetic contribution to the entropy of 
formation for the vacancy in paramagnetic austenite. 



B. Interstitial formation 

All reference states give the (001) dumbbell as the sta- 
ble self-interstitial configuration. The fct afm states sug- 
gest a formation energy of between 3.2 and 3.6 eV, large 
enough to preclude the formation of thermal interstitials 
or Frenkel pairs. Therefore, these defects are only of im- 
portance in irradiated samples. The symmetry breaking 
effect between the [001] and [100] dumbbells is approx- 
imately 0.1 eV i.e. around 3% of the formation energy. 
Experimental evidence in the fee metals aluminium and 
copper indicates the (001) dumbbell as the most stable 
self-interstitial defect™. Other measurements^i^i^ are 
consistent with this conclusion and also suggest the same 
to be true of fee Ni. In summary all experimental evi- 
dence suggests that the (001) dumbbell is the most sta- 
ble self-interstitial in fee metals, with theory in strong 
agreement^iSii^"— here. 

Magnetic moments on the (001) dumbbell atoms were 
severely reduced or even flipped relative to the magnetic 
plane containing the dumbbell in all reference states. Mo- 
ment flips were observed for the [001] dumbbell in the fm- 
HS state and for the [100] dumbbell in the afml state. In 
all other cases the moments were reduced to values in the 
range from 0.1 /is to 0.5 ^b- The influence of the (001) 
dumbbell on the local lattice naturally splits the Inn shell 
into tensile sites, lying within the plane perpendicular to 
the dumbbell axis, and compressive sites lying above and 
below that plane. The effect on the compressive sites 
is by far the most pronounced, exhibiting displacements 
of between 0.2 and 0.3 A consistently for all reference 
states. Magnetic moments on these atoms were consis- 
tently reduced in magnitude by at least 0.33 /is with a 
greatest reduction of 0.92 ^b observed in the afml state 
for a [001] dumbbell. Atoms in the tensile sites at Inn 
relaxed towards the dumbbell centre by between 0.06 and 
0.13 A and exhibited enhanced moments of 0.10 /is in 
the fm-HS state and between 0.21 and 0.44 jib in the afm 
states. 



C. Vacancy migration 

We present estimates of formation energies for the 
transition states involved in vacancy migration and the 
respective barriers heights above the energy of an isolated 
vacancy in TABLE |llll 

The transition states used were those naturally sug- 
gested by symmetry with the migrating atom placed 
half way between the two lattice sites involved. Nudged 
elastic band (NEB) method calculations confirm this to 



Path fee afmD fct afmD fct afml fct fm-HS 

E{ Era Et Era Ef Em E{ Era 



la 


2.717 1.046 2.563 0.743 2.575 0.622 1.826 0.133 


lb 


2.384 0.712 2.867 1.048 N/A (-3.935 -5.627) 


Ic 


2.940 1.268 3.401 1.581 3.677 1.724 N/A 



TABLE III: Formation energies, Et, and estimated migration 
barrier heights, Em, in eV for vacancy point defects in fcc/fct 
Fe. The migration paths are labelled by the two sites involved 
as given in FIG. 2(b) 



be the correct choice for migration within the magnetic 
planes of the afmD state. 

The results for the fct fm-HS state show clear signs of 
instability and are only included here to illustrate this 
point. The broken symmetry of the afm states means 
that there are several non-equivalent barriers. For path 
Ic the initial and final moments on the migrating atom 
are opposite in sign. The constraint of collinear mag- 
netism means the moment must therefore either be zero 
at some point along the path, which is very likely to be 
the transition state in that case given the high energy 
cost of suppressing the moment in Fe, or discontinuously 
flip sign at some point. In our calculations a stable mo- 
ment of zero was ensured by symmetry in the transition 
state, resulting in relatively large energy barriers. How- 
ever, there is no reason to constrain the moment to be 
continuous along the migration path and the inclusion 
of discontinuous flips would very likely lower the barrier 
height and the values in TABLE IIIII should therefore be 
taken as upper limits. 

Along paths la and lb the initial and flnal moments 
for the migrating atom have the same sign and so no mo- 
ment flips are required. In the afml state, however, the 
most stable moment for the migrating atom was found to 
have opposite sign to the initial and final points. A stable 
solution was found with the same sign of moment but this 
was found to be around 0.4 eV higher in energy. Despite 
this complication, our estimate of the barrier height still 
stands if we allow discontinuous moment flips along the 
migration path. Our data shows that moment flips have 
an energy cost of around 0.05 eV and would therefore 
have no effect on barrier height. We therefore estimate 
the relevant barriers for vacancy migration in austenite 
to be in the range from 0.6 to 1.05 eV. Combined with the 
formation energy, this gives good agreement with the ex- 
perimental activation energy for self-diffusion in austen- 
ite which is Qo = 2.945 eV^. 

One final point worth mentioning is that the moments 
on the migrating atoms in the transition state were con- 
sistently enhanced relative to bulk to between 2.51 and 
2.68 hb, unless constrained to be zero. These increases 
exceed those found for atoms Inn to a vacancy defect, as 
would be expected given the larger volume the migrating 
atoms occupy. 



D. Point defect interactions and clustering 

The fate of an irradiated material is initially dependent 
on the interactions between point defects. Calculations 
of such quantities are given in TABLES IIVI andlVl Again, 
the low symmetry of the reference states necessitates cal- 
culation of many configurations, however a clear picture 
emerges from this, that vacancies bind to form di vacan- 
cies with an energy of order 0.1 eV (or in afmD up to 0.2 
eV). This rather weak binding suggests that at elevated 
temperature, divacancies will not be thermodynamically 
stable and nucleation of voids will face a nucleation bar- 
rier. By contrast, interstitials bind strongly into pairs, 
with binding energies of around 1 eV for parallel [001] 
dumbbells on adjacent sites, just as was observed in sim- 
ilar calculations in fee Ni^ where a binding of 0.97 eV 
was found. Such structures can form the nucleus of dis- 
location loops. 

These geometric conclusions can be taken as robust, 
given the good agreement between the two magnetic 
structures. The elastic strain fields are still small (~ 
0.02 eV) for di- vacancies, however, for di-interstitials the 
Eshelby correction lies between -0.2 and -0.3 eV. The ef- 
fect on binding energies is smaller due to a partial cancel- 
lation of correction terms but is still significant, resulting 
in between 0.1 and 0.15 eV increases. However, the sta- 
ble geometry is determined by the differences between 
binding energies, and these converge much more rapidly. 



A-B/Cfg. 



fct afmD 

Ef Eh 



fct afml 

Ef Eh 



Cfg. 



fee afmD 

E{ Eh 



id afmD 

E{ Eh 



fct afml 

Ef Eh 



fct fm-HS 

Ef Eh 



la 3.139 0.205 3.602 0.037 3.843 0.063 3.558 -0.173 

lb 3.288 0.056 3.512 0.127 N/A (-2.586 5.970) 

Ic 3.269 0.075 3.463 0.175 3.860 0.046 N/A 

2a 3.322 0.022 3.702 -0.064 3.883 0.023 (-2.600 5.984) 

2b 3.423 -0.079 3.657 -0.018 3.996 -0.090 3.348 0.036 

TABLE IV: Formation energies, Ef, and, binding energies, 
Eh, in eV for interacting vacancies at up to second near- 
est n eighb our separation. Configurations are labelled as in 
FIG.[2(b)| 



To investigate void formation and dislocation loop nu- 
cleation further we have performed a set of small defect 
cluster calculations. Our choice of vacancy cluster con- 
figurations was motivated by the observation that the 
strongest vacancy-vacancy binding energies are at Inn 
separation (TABLE IIVI) . Our results for clusters con- 
taining up to 6 vacancies are given in TABLE I VII and 
presented graphically along with results for divacancies 
in FIG. 



4(a; 



For dislocation loop nucleation we consider 
clusters of up to 5 [001] dumbbells (axis perpendicular to 
magnetic planes) lying within a single magnetic plane, as 
motivated by the strong Inn binding seen for both afml 
and afmD. The results are presented in TABLE IVIII and 
shown along with the pair binding results in FIG. |4(b)[ 
The first feature which is notable is that both vacancy 
and interstitial defects have a strong tendency to bind 



[001]-[001]/la 5.541 0.850 6.423 0.807 
[001]-[001]/lb 6.584 -0.194 N/A 
[001]-[001]/lc 6.339 0.052 7.301 -0.071 
[001]- [001] /2a 6.556 -0.165 7.789 -0.559 
[001]-[001]/2b 6.398 -0.007 7.344 -0.114 


[001]-[100]/la 5.942 0.569 6.655 0.491 
[001]-[100]/lb 5.629 0.882 N/A 
[001]-[010]/lb unstable N/A 
[001]-[100]/lc 5.867 0.644 6.570 0.576 
[001]-[010]/lc unstable unstable 
[001]-[100]/2a 6.284 0.227 6.987 0.159 
[001]- [010] /2a 6.084 0.428 6.711 0.435 
[001]-[100]/2b 6.101 0.410 7.017 0.129 
[001]-[100]/2c 6.339 0.172 N/A 


[100]-[100]/la 6.417 0.214 7.085 -0.023 
[100]-[010]/la unstable unstable 
[100]-[100]/lb 5.928 0.703 N/A 
[010]-[010]/lb 6.506 0.125 N/A 
[100]-[010]/lb 6.120 0.511 N/A 
[100]-[100]/lc 5.540 1.091 6.539 0.523 
[010]-[010]/lc 6.451 0.181 7.161 0.099 
[100]-[010]/lc 5.988 0.643 6.618 0.444 
[100]- [100] /2a unstable unstable 
[010]-[010]/2a 6.760 -0.129 7.667 -0.605 
[100]- [010] /2a 6.238 0.393 6.913 0.149 
[100]-[100]/2b 6.572 0.059 7.548 -0.486 
[100]-[010]/2b 6.383 0.248 6.846 0.216 



TABLE V: Formation energies, Ef, and, binding energies. Eh, 
in eV for interacting (100) SI dumbbells. Configurations are 
labelled as in FIG. 2(b) with the ordering of specific dumb- 
bells, labelled as A and B, uniquely identifying a configura- 
tion. 



into clusters. As shown in FIG. HI the total binding 
energy increases steadily with the addition of defects. 
This is especially pronounced for the vacancy cluster data 
where, despite the large spread of values for a fixed num- 
ber of vacancies, there appears to be a super-linear in- 
crease in the total binding energy. This implies that, on 
average, the binding of an additional vacancy to an al- 
ready existing cluster is increasing with cluster size. The 
additional binding should eventually tend to the forma- 
tion energy for a vacancy in bulk although it appears 
from the data we are not at this limit yet. 

The most strongly bound tri- vacancy cluster is formed 
by the removal of a tetrahedron of atoms, with an addi- 
tional atom replaced near the centre. There is a strong 
agreement between the afml and afmD states here with 
a total binding energy of 0.58 eV and a magnetic mo- 
ment on the central atom of 3 /ib- Our calculations also 
show that if the central Fe atom is positioned symmet- 



Config. 



fct afmD 

E{ El, 



fct afml 

E[ El, 



(0,5,6) 

(0,9,11) 

(0,1,3) 

(0,1,9) 

(0,5,7) 

(0,1,11) 

(0,5,11) 

(0,3,5) 

(0,1,10) 

(0,5,1) 

(0,5,9) 

(0,5,9,10) + Fe sym 

(0,5,9,10) + Fe asym 

(0,1,2,5) + Fe 



5.388 0.069 
5.094 0.363 
5.179 0.279 
5.134 0.324 
5.348 0.110 
5.112 0.346 
5.264 0.194 
5.307 0.151 
5.128 0.329 
unstable 
unstable 
5.662 -0.204 
4.883 0.575 
4.815 0.642 



5.711 0.148 

5.746 0.113 
N/A 

5.800 0.059 
5.694 0.165 

5.747 0.112 
5.761 0.098 

N/A 
5.761 0.098 

N/A 
5.507 0.352 
6.072 -0.213 
5.283 0.576 

N/A 



(0,5,9,10) 

(0,1,2,5) 

(0,5,6,13) 

(0,9,11,14) 



6.576 0.701 

6.501 0.776 

6.994 0.282 

6.499 0.777 



7.110 0.702 

N/A 
7.406 0.406 
7.496 0.316 



(0,5,6,10,13) 

(0,2,5,6,13) 

(2,5,6,10,13) 

(0,2,5,6,10,13)-^ Fe 



7.866 1.230 

8.220 0.876 

7.507 1.589 

7.404 1.692 



8.337 1.428 

N/A 
8.312 1.451 
8.183 1.582 



(0,2,5,6,10,13) 

(0,1,4,5,7,16) PD 

(0,1,4,5,7,16) SFT 

(0,5,7,10,11,17) SFT 



8.402 2.513 

unstable 
8.573 2.342 
8.751 2.164 



9.200 2.518 
10.532 1.184 
9.366 2.350 

N/A 



TABLE VI: Formation, Ef, and total binding, E^,, energies in 
eV for vacancy clusters. Configurations are identified by the 
lattice sites occupied by vacancies, as numbered in FIG. [3] 
Each section of the table contains configurations with the 
same number of vacancies. Some of the configurations consist 
of vacancy clusters forming symmetrical voids with a single 
Fe atom placed at or near the centre (denoted by +Fe). In the 
case of the tetrahedral (0,5,9,10)-|-Fe cluster stable configura- 
tions were found with the central Fe atom placed symmetri- 
cally (denoted sym) and off-centre along an axis perpendicular 
to the magnetic planes (denoted asym). For the six vacancy 
clusters we distinguish between a (111) planar defect (PD) 
and a stacking fault tetrahedron (SFT) having the same base 
as the planar defect. 



rically, with zero moment, the energy is almost 0.8 eV 
higher, indicating the high cost of suppressing the Fe 
moment. Such a large value for the moment is consistent 
with earlier observations that the local moment increases 
along with the local volume occupied by an atom (unless 
constrained to zero by symmetry) , converging asymptot- 
ically to the moment on the free atom. Calculations in 
bcc Fe are in strong agreement with this finding^^. 

Significant differences between the reference states 
were observed for the binding energies of the less well 
bound tri-vacancy clusters. A particular case worthy 




FIG. 3: Lattice site numbering used to identify clusters of 
point defects and solutes. The afmD magnetic state is shown 
to identify its clusters unambiguously but the numbering is 
equally valid for the other magnetic states considered here 
(where no ambiguity exists). 

Config. fct afmD fct afml 

Ef Eh Ef Eh 



(0,5,6) 


8.048 


1.538 


9.580 


1.265 


(0,5,7) 


7.897 


1.689 


9.168 


1.677 


(0,5,6,7) 


10.583 


2.198 


13.006 


1.454 


(0,5,6,13) 


9.820 


2.962 


11.867 


2.593 


(5,6,7,8) 


13.055 


-0.274 


15.628 


-1.168 


(0,5,6,7,8) 


13.061 


2.916 


16.651 


1.424 



TABLE VII: Formation, Ef, and total binding. Eh, energies 
in eV for clusters of [001] SI dumbbells. Configurations are 
identified by the lattice sites occupied by dumbbells, as num- 
bered in FIG. El 



of mention is configuration (0, 5, 9) (see TABLE IVII 
and FIG. [3]), which is stable in the afml state but found to 
be unstable for afmD, where it relaxes to one of the tetra- 
hedral arrangements. Despite this, however, the next 
most strongly bound tri-vacancy cluster in both cases has 
a binding energy of approximately 0.35 eV, clearly indi- 
cating the additional stability of the tetrahedral arrange- 
ment. It is well established that stacking fault tetrahedra 
(SFT) are stable vacancy configurations in fee and this 
cluster could be regarded as the smallest possible such 
object. 

Another vacancy can be absorbed to form a tetrahe- 
dron, with an additional binding of approximately 0.1 eV 
but this is significantly lower than the binding energy per 
defect of the most stable tri-vacancy so this may prove a 
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FIG. 4: Total binding energies, Eh, in eV for (a) vacancy 
clusters and (b) self-interstitial dumbbell clusters. Data has 
been shifted horizontally for clarity. Configurations of defects 
are considered clusters when all defects can be connected by 
chains of Inn bonds. 



bottleneck against forming three dimensional voids. 

There is, however, a rapid increase in the stability of 
vacancy clusters above this point. The addition of a fifth 
vacancy increases the binding energy by around 0.7 eV 
on average. The most stable five vacancy configuration 
considered here is formed by adding a single Fe atom to 
the centre of an octahedral vacancy cluster. The local 
arrangement of first neighbours to this central atom is 
bcc-like. There is a significant relaxation of these first 
neighbours towards the central atom but the effective bcc 
lattice parameter is still greater than that for bcc Fe at 
equilibrium and significant tetragonal distortion remains. 
The moment on the central atom was found to align with 
its neighbours in the afml state i.e. locally ferromagnetic 
with moments between 2.3 fiB and 2.4 fj,B- In afmD 



similarly enhanced moments relative to the bulk were 
found but the afmD magnetic order remained the stable 
arrangement. 



The jump to six vacancy clusters is again associated 
with a large increase in total binding energy. Of partic- 
ular interest is the extremely stable octahedral arrange- 
ment of six vacancies with a total binding energy of 2.5 
eV in both the afml and afmD states, which is 0.8 eV 
higher than the binding energy of any five vacancy clus- 
ter. The six vacancy SFT configurations show similarly 
high binding energies. These configurations are the next 
size of SFT up from the elementary example discussed 
earlier and can be formed by removing six atoms in a Inn 
triangular arrangement from a (111) plane in fee with a 
subsequent and significant relaxation of the four atoms 
forming a tetrahedron directly above this plane down- 
ward to fill the void. This type of defect represents an 
alternative to the formation of three dimensional voids, 
best represented in our calculations by the octahedral 
cluster. We cannot clearly distinguish between these two 
on energetic grounds. One interesting aside here is that 
we found the six vacancy (111) planar defect to be meta- 
stable in the afml state whereas this arrangement relaxed 
to the SFT in afmD. This mirrors what was found for 
the three vacancy (111) planar (0,5,9) cluster. In both 
of these cases, however, the SFT was found to be signif- 
icantly more stable. 



The Eshelby corrections to the formation energy of 
most of our vacancy cluster calculations are only a few 
hundredths of eV, with some notable exceptions. Cor- 
rections for five vacancy clusters are at most -0.13 eV, 
for the octahedral cluster the correction is -0.06eV and 
for the stacking fault tetrahedra the correction is -0.2eV. 
These corrections do not affect our conclusions in any 
significant way. 



For interstitial clusters we also find a strong tendency 
to bind and align in parallel to form (001) proto-loops 
(c.f FIG. 4(b) ), again in agreement with ab initio results 
in fee Ni^ where a tri-interstitial binding energy of 1.71 



eV was observed. The picture is complicated by repul- 
sion between 2nn dumbbells, which is particularly strong 
for the fct afml state. The most strongly bound planar 
clusters are therefore those which maximise the ratio of 
Inn to 2nn bonds. The precise energetics and geometry 
of these results should be treated with care, since our re- 
sults sample a few of the very large number of possible 
configurations. Furthermore, the Eshelby correction to 
the binding energy was found to be around -leV for the 
largest cluster (as confirmed by a direct constant pres- 
sure, P = GPa, calculation), a value that is comparable 
with the binding itself. Even so, this correction can only 
increase the binding energy, enhancing the driving force 
for clustering. 



E. Pair bond models for defect clusters 

We have attempted to model the binding energy of 
pairs of defects and defect clusters with a linear pairwise 
bonding model up to second nearest neighbours, 



4'"^ (^1,^2) 



PlTli +P2n2, 



(3) 



where ni and ^2 are the total number of Inn and 2nn 
bonds respectively between defects in the cluster and 
Pi and p2 are the corresponding fit parameters. This 
model does not distinguish between the distinct neigh- 
bour bonds resulting from the symmetry breaking effects 
of the magnetic states and tetragonal distortion and fits 
using this model therefore provide an effective averaging 
over distinct bonds. Including these symmetry breaking 
effects was found to significantly improve the agreement 
between model and data but at the expense of consider- 
ably more parameters. 

The nature of the model also means that it is only di- 
rectly applicable to clusters where the individual defects 
can be assigned to single lattice sites. This is not the 
case for the tetrahedral and octahedral voids with a sin- 
gle Fe atom at the centre and for the six vacancy stacking 
fault tetrahedra and these configurations have not been 
included in the modelling. Fits were performed both to 
the afml and afmD data sets individually and to their 
combined data, averaging over the two sets in the pro- 
cess. For vacancy clusters the inclusion of 2nn bonds in 
the model made little difference to the results and have 
therefore been omitted (effectively setting p2 =0). The 
resulting fit parameters are given in TABLE [Villi and the 
corresponding model values are compared with the data 
in FIG. El 
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Parameter afmD afml 



both 



Vacancy cluster results 



pi/eV 0.175 0.141 



0.162 



Self-interstitial dumbbell cluster results 
pi/eV 0.841 0.829 0.841 

P2/eV -0.142 -0.441 -0.275 

TABLE VIII: Fit parameters, Pi, for the pair bond model 
from fits to defect cluster binding energies for single magnetic 
data sets or combined fits to both data sets as indicated. 

The fit parameters for these models have the units of 
energy and can be thought of as effective or averaged 
binding energies for pairs of defects. Comparison with 
the pair binding energies shows that this is generally true, 
particularly for interstitials. The afml vacancy cluster 
parameter, pi = 0.141 eV, does, however, overestimate 
the pair binding but should be thought of as an effective 
value which best represents the interactions in all clusters 
included in the fit. 



From FIGURES 5(a) and 5(b) it is clear that there is a 



significant horizontal spread in the vacancy cluster ab ini- 
tio data for a fixed number of bonds and therefore model 




• fct afmD single 
■ fct afml single 
O fct afmD both 
□ fct afml both 

„ (m) (d) 



Eb /eV 

(c) Dumbbell cluster 



FIG. 5: Plots of binding energies from the pairwise bond 



r(m) 



p(d) 



model, E^^ , versus ab initio data, E^ , for defect clusters in 
the fct afml and afmD states. The results of fits performed to 
data from a single magnetic state or simultaneously to both 
are shown, as labelled. The simultaneous fit results have been 
separated into those for afml and afmD states in order to al- 
low the effect of a combined fit to be directly compared with 



p(m) 



p(d) 



the single fit results. The E^ — E^ line is included to 
indicate a perfect fit. Horizontal lines indicate representative 
error bars for the averaging over both magnetic states in the 
combined fits. 
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value, just as there was for a fixed number of vacancies 
in FIG. 4(a) Tfiis is a direct result of symmetry breaking 



effects in the reference states, which the model does not 
incorporate and effectively averages over. There is, how- 
ever, a broadly linear trend in the data that the model 
is able to capture. This is even true of the combined fit 
to both data sets where the model values lie generally 
within the error bars resulting from averaging over dis- 
tinct bonds and magnetic states. It is worth mentioning 
that for some of the data points the model values for the 
combined fit were not significantly different from the sin- 
gle fits and these points in the plots are obscured as a 
result. Despite the general agreement of the models they 
do tend to overestimate the average binding of smaller 
clusters and underestimate the binding of the largest (oc- 
tahedral) cluster. This is certainly attributable to the 
exceedingly strong binding of the octahedral configura- 
tion which may well define a limit on the applicability of 
this model since still larger clusters are likely to be even 
more strongly bound. 

The agreement between models and data for intersti- 
tial clusters (shown in FIG. 5(c)| is significantly better 
than for vacancies. This is partly due to our consider- 
ation of only planar clusters but also demonstrates the 
existence of clear trends in the data. The primary differ- 
ence between the two magnetic states is the significantly 
larger repulsion between dumbbells at 2nn separation for 
the afml state. Data for clusters where 2nn interactions 
are significant therefore shows significant spread. Model 
values for the combined fit offer an effective averaging 
over these differences and provide our best guess for a 
predictive model of planar defects in austenite. 



F. Comparison with bcc fm Fe 

As a final means to summarise our findings in pure 
Fe we make a direct comparison of our data for the afm 
states of austenitic Fe to the ferritic (bcc fm) ground 
state, as shown in FIG. [51 

We have included data for the formation of single va- 
cancy and interstitial defectsi^, the di-vacancy binding 
energy and the binding of the most stable (dumbbell) 
interstitial defect into planar clusters^. Each quantity 
plotted is uniquely defined in the bcc fm state with the 
vertical spread in fee values coming from the symmetry 
breaking effects in the magnetic states considered here. 

The data clearly shows a high level of similarity be- 
tween the bcc and fee results at the energy resolution of 
the plot. Each quantity occupies a clearly defined energy 
range for the corresponding process e.g. di-vacancy bind- 
ing from 0.0 to 0.2 eV and di-interstitial binding from 
0.8 to 0.9 eV. It also demonstrates the generally good 
level of agreement between the afml and afmD states. 
One noticeable difference is the generally lower forma- 
tion energies for point defects in austenite, particularly 
for afmD. Despite this, the binding energies show a gen- 
erally good agreement. Overall this raises the question 




^/ Interstitial Formation 

Quad-Interstitial Binding 
Q Vacancy Formation 
□ Tri-Interstitial Binding 
Di-Interstital Binding 

Di-Vacancy Binding 

12 3 4 

E(bcc) / eV 



FIG. 6: Comparison of point defect formation and binding 
energies between austenitic and ferritic Fe. Interstitial for- 
mation energies for the two datasets have been associated 
with one another in order from lowest to highest energy. The 
black line is included to indicate an exact agreement. 



of whether such a crystal structure independence exists 
generally in metals or is particular to Fe. 



V. DILUTE NI,CR RESULTS 

Gommercial austenitic stainless steels typically contain 
Cr and Ni as major alloying elements. As a first step 
we have examined the behaviour of these atoms in the 
dilute limit. Our results for single substitutional solutes 
are given for the four reference states considered here 
in TABLE HXl 

Config. fee afmD fct afmD fct afml fct fm-HS 

Ei fi E[ fi Ef ^ Ei fi 

Sub. Ni -0.033 -0.08 0.084 0.04 0.167 -0.29 -0.053 0.70 

Sub. Cr 0.106 0.53 0.268 0.85 0.047 1.07 0.036 -2.17 

TABLE IX: Formation energies, Ef, in eV and magnetic 
moments, fi, in fiB for substitutional Ni and Cr solutes in 
austenitic Fe. The sign of the moments indicates whether 
there is alignment (positive) or anti-alignment (negative) with 
the moments of the atoms in the same magnetic plane. 

The substitutional formation energy of both Ni and Gr 
in austenitic Fe is small but positive in the fct afm ref- 
erence states indicating a weak tendency for phase seg- 
regation, but only at temperatures way below the sta- 
bility limit of austenite. In the fm-HS state we find Ni 
to be soluble whereas Cr is not and whilst the energies 
involved are small this trend is opposite to that seen in 
similar calculations in bcc fm Fe^. The magnetic mo- 
ments on the solute atoms exhibit the usual behaviour 
of Ni to be ferromagnetic to its neighbours and Gr to 
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be anti-ferromagnetic. In the afmD state the Ni mo- 
ment is heavily suppressed relative to its pure reference 
state value (/i^f — 0.59 [Ib) whereas in afml the mo- 
ment is aligned with the majority of its Inn Fe atoms 
although still reduced in magnitude. The moment on a 
single Cr atom in the afm states remains comparable to 
the reference state value (A^cr ~ 0.89 /is) and even shows 
enhancement in the afml state. Enhanced moments are 
also seen for both Ni and Cr solutes in the fm-HS state, 
exhibiting alignment and anti-alignment with Fe respec- 
tively, just as was observed in bcc fm Fe^. 

Bond lengths from single Ni and Cr solutes to their 
nearest neighbour shells differ from pure Fe by at most 
0.05 A. This 1 to 2 % effect, exists only for the Inn shell 
and decays rapidly with distance. A study of the magni- 
tude and directions of bond length changes is complicated 
by the symmetry breaking effects in the reference states. 
We consider instead the changes in the lattice parameters 
of the unit cell surrounding a single substitutional solute, 
which naturally distinguishes between effects within and 
perpendicular to the magnetic planes. The build up of 
stresses on the supercells were found to be consistent with 
the changes in the lattice parameters. The influence of Ni 
and Cr solutes was found to be very similar within the 
afm reference states but showed differences in the fm- 
HS state. For the fct afm states both solutes increased 
the lattice parameter within a magnetic plane by 0.04 A 
in afmD and 0.05 A in afml. The influence on the out 
of plane lattice parameter distinguished between afmD, 
which showed a contraction of 0.03 A, and afml, where 
an almost negligible increase was found for Ni and a small 
contraction for Cr of 0.016 A. The only significant change 
in the fm-HS state was a contraction of 0.04 A in the out 
of plane lattice parameter for Ni. 

The influence of Ni and Cr solutes on the magnetic 
moments of surrounding Fe atoms is more pronounced. 
For the Inn shell changes of 0.1 /is are typical but were 
found to be as high as 0.16 /is for Ni in the afml state. 
No significant (> 1 %) moment changes were observed 
for the fm-HS state at higher orders. In the fct afmD 
state moments 2nn to Ni showed perturbations of up to 
0.06 yuB and were negligible above that. The influence of 
Cr was similar in magnitude at 2nn but persisted out to 
the 4nn shell. In the fct afml state both Ni and Cr solutes 
exerted a more pronounced influence out to the 4nn shell 
than in other states. For Ni, moments on atoms at 2nn 
differed by up to 0.12 \xb (i.e. 8 %) from the pure Fe value 
and differences of up to 0.05 [Ib were observed at 4nn. 
For Cr, moment changes in the 4nn shell were the largest 
measured for that separation at 0.07 \JbB^ exceeding those 
at Inn separation, which were up to 0.05 /is. 

Overall we conclude that a significant contribution to 
the interactions of Ni and Cr solutes with defects and 
other solutes in austenitic Fe, especially at longer range, 
will come from their magnetic interactions. Volume- 
elastic contributions will be smaller and should be similar 
for Ni and Cr in the fct afmD and fct afml states. 

Interactions between pairs of Ni and Cr solutes, given 



Config. 
A-B/cfg 



fee afmD 

El Eh 



fct afmD 

Ei Eh 



fct afml 

Ei Eh 



fct fm-HS 

Ef Eh 



Ni-Ni/la 
Ni-Ni/lb 
Ni-Ni/lc 
Ni-Ni/2a 
Ni-Ni/2b 



-0.089 
-0.061 
-0.022 
-0.101 
-0.020 



0.024 
-0.004 
-0.043 
0.036 
-0.045 



0.113 
0.141 
0.182 
0.151 
0.182 



0.056 
0.027 
-0.014 
0.018 
-0.014 



0.204 0.110 

N/A 
0.280 0.034 
0.267 0.067 
0.257 0.058 



-0.124 0.018 
-0.126 0.020 

N/A 
(-5.853 5.748) 
-0.086 -0.020 



Cr-Cr/la 
Cr-Cr/lb 
Cr-Cr/lc 
Cr-Cr/2a 
Cr-Cr/2b 



0.257 
0.218 
0.305 
0.170 
0.214 



-0.044 
-0.005 
-0.093 
0.043 
-0.001 



0.562 
0.547 
0.633 
0.512 
0.546 



-0.027 
-0.012 
-0.098 
0.023 
-0.011 



0.173 -0.062 

N/A 
0.183 -0.071 
0.086 0.008 
0.121 -0.010 



0.313 -0.241 
(-5.742 5.814) 

N/A 
0.102 -0.030 
0.077 -0.005 



Ni-Cr/la 
Ni-Cr/lb 
Ni-Cr/lc 
Ni-Cr/2a 
Ni-Cr/2b 
Ni-Cr/2c 



0.090 
0.092 
0.095 
0.095 
0.095 
0.098 



-0.016 
-0.018 
-0.021 
-0.021 
-0.021 
-0.024 



0.327 
0.363 
0.325 
0.367 
0.367 
0.359 



0.025 
-0.012 
0.027 
-0.016 
-0.016 
-0.008 



0.162 0.051 

N/A 
0.215 -0.002 
0.236 -0.022 
0.239 -0.025 

N/A 



0.012 -0.028 
(-5.750 5.734) 

N/A 

-0.029 0.012 

0.013 -0.029 

N/A 



TABLE X: Formation energies, E{, and, binding energies, 
Eh, in cV for substitutional Ni and Cr solutes in Fe. Solute 
pair co nfigura tions are labelled according to the definitions 
in FIG. 



2(b) 



The afml data has been produced with a 3 
Monkhorst-Pack kpt sampling, except for Ni-Cr 2b. Negative 
formation energies for the fct fm-HS state have once again 
been included (in brackets) only to illustrate the instability 
of the reference state. 



in TABLE |Xl are weak in general. The only non- 
negligible attraction observed here is between pairs of 
Ni atoms of up to 0.1 eV at Inn separation and at 
around 0.06 eV for the afml state at 2nn. This ex- 
tended attraction may well result from the particularly 
strong and long ranging infiuence of Ni on neighbour- 
ing magnetic moments in that state. Interactions be- 
tween pairs of Cr atoms are repulsive at Inn separation 
and found to be particularly strong for the fm-HS state 
where E^y — —0.24 eV. A very similar negative binding 
was found for Cr-Cr pairs in bcc fm Fe^ at Inn separa- 
tion. The binding of Ni-Ni pairs was found to be almost 
negligible in bcc Fe, in contrast to the modest attraction 
seen in afm states here, but consistent with our findings 
for the fm-HS state. Interactions between pairs of Ni and 
Cr are mostly negligible save for some signs of an attrac- 
tion at Inn in the afm states. Based on these solute pair 
interactions we suggest that the amount of short range 
order in alloys at typical operating temperatures for nu- 
clear applications will be small but with a tendency for 
Ni rich clustering and Cr-Cr repulsion. 

We have performed a small set of solute cluster calcu- 
lations, containing up to 5 solute atoms, in order to fur- 
ther investigate the trends seen for pairs of solutes and 
to investigate whether a simple pair interaction model is 
consistent with the data. The results of our calculations, 
which were performed only for the fct afm states, are 
given in TABLE |Xll 
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config. 




fct afmD 


fct afml 














Ei 


Sb 


Et 


£b 


> 








Ni 


(0,5,7) 
(0,5,6) 
(0,5,7) 
(0,5,6) 




0.141 
0.124 
0.795 
0.815 
0.549 


0.111 
0.128 
0.008 
-0.012 
0.071 


0.272 
0.231 
0.266 
0.256 
0.155 


0.227 
0.269 
-0.125 
-0.115 
0.105 


0.5 
0.4 
0.3 

OT 


' 1 ' 1 


Ni 
Cr 
Cr 


• 
- ■ 


fct afmD 
fct afml 


Ni:(0), Cr:(5,6) 






Cr:(0), Ni:(5,6) 




0.361 


0.075 


0.202 


0.179 


E 








Ni;(5,6,14,15) 




0.336 


-0.000 


0.430 


0.236 


V 


U.l 


~ 


■ 


Cr:(5,6,14,15) 




1.085 


-0.015 


0.251 


-0.063 







- 


•sit 


Ni:(0), Ni:(5,6,14, 


15) 


0.326 


0.094 


0.367 


0.466 




-0 1 




J^ 


Ni:(0), Cr:(5,6,14, 


15) 


1.138 


0.016 


0.357 


-0.003 






m^ 


/V 


Cr:(0), Ni:(5,6,14, 


15) 


0.559 


0.045 


0.422 


0.291 




-0.2 


-^ 


III, 




OT 


-0.1 


Cr:(0), Cr:(5,6,14, 


15) 


1.386 


-0.048 


0.422 


-0.187 











0.1 

. (d) 



0.2 

/eV 



0.3 



0.4 



0.5 



TABLE XI: Formation, .Bf, and binding, _Eb, energies in eV 
for solute clusters. Configurations are identified by listing 
the lattice sites occupied by Ni and Cr solutes, as numbered 
in FIG. El 



It is immediately apparent that the most strongly 
bound solute clusters contain predominantly Ni atoms. 
More particularly these clusters contain a majority of Ni- 
Ni Inn and 2nn pair bonds and Ni-Cr Inn bonds indicat- 
ing these results are consistent with the binding energies 
of pairs of interacting solutes and that a pair interaction 
model is justified. In order to quantify this claim we have 
performed fits to our solute cluster and solute pair bind- 
ing energies using a pair bond model. The form of this 
model is similar to that used for defect clusters but now 
we have to count the numbers of Ni-Ni, Ni-Cr and Cr-Cr 
bonds independently, resulting in a 6 parameter model: 



E. 



(m) _ Ni-Ni^Ni-Ni , Ni-Ni^Ni-Ni 



Vi 



V2 



Pi 



Cr-Ci„ Cr-Cr 



n. 



P2 



Cr-Cr^ Cr-Cr 



Tin 



Pi 



Ni-Cr^Ni-Cr 



n. 



-P2 



Ni-Cr^Ni-Cr 



Tin 



(4) 



Fits to the fct afmD dataset with this model showed 
good agreement within errors. For the afml dataset the 
agreement was also generally good but there were clear 
outliers from the fit. Further inspection showed that 
these outliers came exclusively from configurations con- 
taining Ni atoms in the (5,6,14,15) positions i.e. those 
with Ni-Ni 3nn bonds. A calculation of the relevant Ni- 
Ni 3nn binding energy in afml gave a value of 0.055 eV. 
Within a pair bond model this result is consistent with 
the relatively high binding energies seen in the outliers 
and we have therefore included an extra term to include 
contributions from Ni-Ni 3nn bonds in the afml fit. We 
have not attempted a combined fit to the afml and afmD 
datasets because the magnitude and range of the solute- 
solute interactions is too distinct and the number of fit 
parameters differs between the sets. The results for the 
6 parameter fct afmD and 7 parameter fct afml fits are 
compared with the data in FIG. [7] and the fit parameters 



(a) Solute cluster 




E, /eV 

b 
(b) Solute cluster zoomed 



FIG. 7: Plots of binding energies from the pairwise bond 



model, E, 



(m) 



versus ab initio data, E, 



(d) 



for solute clusters 



in the fct afml and afmD states, 
included to indicate a perfect fit. 



The E, 



(m) 



= eI'^^ line is 



given in TABLE En 



7(a) 



As was noted earlier the data presented in FIG 
shows a generally good agreement between model and 
data for solute clusters. The seemingly greatest disagree- 
ment is present in the central section of the graph, as 
shown in FIG. 7(b) which primarily contains binding en- 
ergies for pairs of solutes. All of the significant outlying 
points, however, come from symmetry breaking effects in 
the reference states, resulting in a spread of data values 
corresponding to a single model value, just as was seen 
for the vacancy cluster data. Our models agree with this 
data in the sense that the model value for a particular 
configuration lies within the spread of data, effectively 
finding an average value that is most consistent with all 
of the cluster data. The generally good agreement be- 
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Parameter afmD afml 



pT 


-^VeV 


0.046 


0.096 


p? 


-°7eV 


-0.014 


-0.058 


pT 


-^VeV 


0.030 


0.055 


pf 


-^VeV 


0.005 


0.047 


pf 


-^VeV 


-0.003 


-0.022 


pT 


-^VeV 


-0.014 


-0.031 


pf 


-^'/eV 


— 


0.043 



Solute/cfg fee afmD fet afmD fet afml 

Et Em E{ Em E{ Em 



TABLE XII: Fit parameters, p- ~ , for the pair bond model 
from fits to solute cluster binding energies for the afml and 
afmD states. 



tween solute pair and cluster data strengthens our sug- 
gestion that the amount of short range ordering in alloys 
will be weak with some tendency to form Ni rich clusters. 
Our results do not, however, rule out the possibility of 
complex many body effects in concentrated alloys. 



A. Defect-Solute interaction 

We present our results for the binding of Ni and Cr 
solutes to a single vacancy defect in TABLE [XIIII and for 
their vacancy mediated migration in TABLE IXIVl 



A-B/efg 



fee afmD 

El E-t, 



fet afmD 

El -Eb 



fet afml 

El E-b 



V-Ni/la 
V-Ni/lb 
V-Ni/lc 

V-Ni/2a 
V-Ni/2b 

V-Ni/2c 



1.596 
1.626 

1.675 
1.614 
1.687 
1.691 



0.043 
0.013 
-0.036 
0.026 
-0.048 
-0.051 



1.847 
1.876 
1.887 
1.906 
1.914 
1.909 



0.056 
0.027 
0.016 
-0.003 
-0.011 
-0.005 



2.031 0.089 

N/A 
2.078 0.042 
2.109 0.011 
2.130 -0.010 

N/A 



V-Cr/la 
V-Cr/lb 
V-Cr/lc 
V-Cr/2a 
V-Cr/2b 
V-Cr/2c 



1.866 
1.845 
1.866 
1.803 
1.818 
1.814 



■0.087 
■0.066 
■0.088 
■0.025 
■0.039 
■0.036 



2.083 
2.161 
2.177 
2.103 
2.153 
2.091 



0.004 
-0.075 
-0.091 
-0.017 
-0.066 
-0.004 



1.970 0.030 

N/A 
2.079 -0.079 
2.052 -0.052 
2.075 -0.075 

N/A 



TABLE XIII: Formation energies, Ef, and, binding energies, 
-Bb, (in eV) for a vacancy defect to substitutional Ni and C r 
solutes in Fe. Configuration labelling defined in FIG. 2(b) 



We find that Ni binds to a vacancy but by no more 
than 0.1 eV at Inn and shows no sign of interaction at 
2nn. In contrast the Cr-vacancy interaction is repulsive 
overall, even at 2nn, with (negative) binding energies as 
low as -0.091 eV. Calculations in bcc Fe^®, in contrast, 
show that both Cr and Ni bind to a vacancy: Cr by 0.2 
eV at 2nn separation and Ni by 0.07 eV at Inn. The 
lack of any strong tendency for vacancy binding suggests 
that the rate of microstructural evolution and of creep 
should be relatively unaffected, at least in dilute alloys 



Ni/la 2.901 1.305 2.738 0.891 


3.014 0.983 


Ni/lb 2.509 0.883 3.049 1.172 


N/A 


Ni/le sym 2.608 0.932 3.066 1.179 


3.441 1.363 


Cr/la 2.712 0.846 2.643 0.560 


2.705 0.735 


Cr/lb 2.417 0.572 2.903 0.742 


N/A 


Cr/lc sym 2.523 0.657 3.021 0.844 


3.101 1.022 



TABLE XIV: Formation energies, Ei, and estimated migra- 
tion barrier heights, Em, in eV for vacancy mediated solute 
migration. The transition state configurations for the migra- 
tion are taken to be when the solute atom is midway between 
the two lattice sites involved in the migration, as labelled 
in FIG. [2(b)] 



under irradiation. There is, however, strong experimen- 
tal evidence that increasing Ni content suppresses void 
formationSi^. The inclusion of small quantities (< 1 at. 
%) of oversized solutes, such as Zr and Hf, in austenitic 
Fe-Cr-Ni alloys were also found to significantly suppress 
void formation and radiation-induced segregation (RIS) 
at grain boundaries^"—. Kato^^iSL suggested the posi- 
tive binding of vacancies to the relatively immobile over- 
sized solutes as a mechanism by enhancing recombina- 
tion and inhibiting vacancy diffusion, which is supported 
by the modelling of Stepanov^^ in the case of RIS. The 
vacancy-Ni binding observed in this work is small but is 
likely to be cumulative and should therefore not be over- 
looked as a contributory mechanism for void suppression 
at higher Ni concentrations. 

The vacancy-mediated migration energies (in TA- 
BLE IXIVp show a distinct difference between Ni and 
Cr. While the lowest Ni migration barrier is between 
0.9 and 1.0 eV, that for Cr is between 0.6 and 0.7 eV. 
These values span the migration barrier heights for Fe 
self-diffusion (TABLE UTTI for all but the Ic path, pre- 
sumably because of the higher cost of suppressing the 
Fe moment to zero, as discussed earlier. The very same 
ordering of migration barrier heights was found for Cr, 
Fe and self-diffusion in fee Ni in ab initio studies by Do- 
main and Becquart^i^ and by Tucker— i^ with barrier 
heights of 0.8, 0.95 and 1.05 eV for Cr, Fe and Ni re- 
spectively. Just such a preferential diffusion is posited 
as the primary mechanism for RIS in Fe-Cr-Ni austenitic 
alloys by the inverse Kirkendall effect^'^--^™, leading to Cr 
depletion and Ni enhancement at grain boundaries and 
other defect sinks. 

We present our results for interstitial Ni and Cr so- 
lutes in TABLE IXVl It is clear from the data that Ni 
is strongly repelled from interstitial sites, negative bind- 
ing showing a direct preference for self-interstitial defects 
and substitutional Ni. A small but positive binding was 
found for the FeNi [001] dumbbell but only for the fet 
afmD state. While interesting, this result is at odds with 
the other reference states so cannot be taken as a general 
conclusion for austenite. The generally repulsive trend 
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Defect 



fee afmD 

E[ El, 



fet afmD 

El El, 



id afml 

El El, 



Oeta Ni (1) 
Tetra Ni uu (2) 
Tetra Ni ud (3) 
[110] Ni erow. (4) 
[Oil] Ni erow. uu (5) 
[Oil] Ni erow. ud (6) 
[100] FeNi dumbbell 
[001] FeNi dumbbell 
[001] NiFe dumbbell 
[111] FeNi dumbbell 
[111] NiFe dumbbell 



unstable 
4.314 -0.765 
3.766 -0.467 

unstable 
4.116 -0.378 
4.006 -0.164 
3.416 -0.471 
3.069 -0.311 
3.097 -0.340 
3.431 -0.132 
3.701 -0.153 



unstable 
4.447 -0.499 
4.086 -0.338 

unstable 
4.507 -0.167 
4.296 -0.044 
3.717 -0.317 
3.215 0.065 
3.469 -0.190 
3.740 0.007 
3.989 -0.041 



4.920 -0.400 

N/A 
5.180 -0.691 

N/A 
5.336 -0.370 
4.112 -0.414 
4.116 -0.334 

as FeNi 
4.378 0.111 

as FeNi 



Octa Cr (1) 
Tetra Cr uu (2) 
Tetra Cr ud (3) 



unstable 
3.550 0.137 
3.305 0.133 



unstable 
3.785 0.347 
3.580 0.351 



4.177 0.223 

N/A 
4.305 0.064 



[110 


Cr crow. (4) 


unstable 


unstable 




[Oil 


Cr crow, uu (5) 


3.740 


0.137 


4.166 0.357 


N/A 


[Oil 


Cr crow, ud (6) 


3.691 


0.290 


4.025 0.411 


4.606 0.259 


[100 


FeCr dumbbell 


3.050 


0.034 


3.385 0.198 


3.583 -0.005 


[001 


FeCr dumbbell 


2.933 


-0.036 


3.184 0.279 


3.592 0.070 


[001 


CrFe dumbbell 


2.850 


0.047 


3.267 0.196 


as FeCr 


[111 


FeCr dumbbell 


3.414 


— 


3.739 — 


4.387 0.219 


[111 


CrFe dumbbell 


unstable 


unstable 


as FeCr 


[100 


NiNi dumbbell 






3.897 -0.413 4.620 -0.755 


[100 


NiCr dumbbell 






3.643 0.024 


3.939 -0.194 


[100 


CrCr dumbbell 






3.826 0.025 


3.930 -0.305 


[001 


NiNi dumbbell 






3.438 -0.075 


4.160 -0.210 


[001 


NiCr dumbbell 






3.344 0.203 


3.943 -0.114 


[001 


CrNi dumbbell 






3.162 0.385 


as NiCr 


[001 


CrCr dumbbell 






3.536 0.195 


3.931 -0.269 



TABLE XV: Formation energies, Ei, and binding energies, 
-Eb, in eV for interstitial defects containing Ni and Cr solute 
atoms. Binding is defined relative to the self-interstitial de- 
fect of the same type and substitutional solute atoms in a 
non-interacti ng st ate. The numbering of interstitial defects 
is as in FIG. 2(a) Dumbbell configurations are identified by 



their direction axis and by their composition, in order along 
that axis, with their centres considered to lie on the lattice site 



identified in FIG. 2(a) All (110) mixed dumbbells were found 
to be unstable. In all of the [111] FeNi mixed dumbbell cal- 
culations the Fe atom relaxed to the nearest tetrahedral site 
and the Ni to a substitutional site at Inn to the Fe. Bind- 
ing energies have been calculated accordingly. Calculations of 
configurations shown to be unstable for afmD have not been 
performed for afml and the column entries left blank. Mixed 
dumbbell configurations have not be calculated for fee afmD. 



is also supported by observations that in mixed dumb- 
bell configurations the Ni atom is generally closer to the 
lattice site than Fe. For the [111] FeNi and [111] NiFe 
dumbbells this asymmetry is so pronounced that the con- 
figuration must be considered as a tetrahedral Fe inter- 
stitial with a substitutional Ni atom at Inn separation 



and the binding energies given in the table have been 
calculated accordingly. Although relatively stable, these 
states are always less stable than mixed (100) dumbbells 
and are important only as possible intermediate states in 
the rotation, migration and disassociation of the stable 
interstitials. 

In contrast, Cr generally shows positive binding to in- 
terstitial sites, an effect that is particularly prevalent in 
the afmD states. In mixed dumbbell configurations the 
trend is for Cr atoms to be farther from the lattice site 
than Fe. This was particularly pronounced for the [111] 
CrFe dumbbell in the afmD state which relaxed to the 
tetrahedral Cr ud configuration. Once again the mixed 
(100) dumbbell configurations are the most stable. 

Our calculations for doubly mixed (100) dumbbells, 
show a clear distinction between the fct afml and 
afmD states. We find that NiCr dumbbells exhibit the 
strongest binding in afmD and while the same configu- 
rations are the most stable for afml the interaction is 
repulsive. One result that the two states do have in com- 
mon, however, is that in NiCr dumbbells the Ni atom is 
generally closest to the lattice site and the effect is more 
pronounced that for mixed FeNi dumbbells. 

The emerging picture here is of a general order of pref- 
erence for the different atomic species to be found within 
overcoordinated defects with Cr being the most stable, 
followed by Fe and finally Ni. There is evidence of the 
same ordering in ab initio studies of bcc fm Fe where Cr 
shows significant binding to overcoordinated deiectiM^ 
and there is a repulsive interaction between Ni and the 
most stable dumbbell^^. Ab initio studies of dilute Cr 
and Fe in fee Ni^'^'^'^ also show that Cr binds strongly to 
the (100) dumbbell in mixed and doubly mixed forms but 
Fe shows little interaction. This preferential binding of 
Cr to overcoordinated defects should result in a positive 
association of Cr with the interstitial defect fiux in irra- 
diated environments^i^'^. This is an interesting result 
as it is in the opposite sense to the flux resulting from 
the vacancy mechanism of the inverse Kirkendall effect. 

The same interstitial mechanism was proposed for Ni 
as a possible explanation for RIS in Fe-Cr-Ni austenitic 
alloys, e.g. the work of Watanabe''^ where a binding of 
0.75 eV was suggested. The inclusion of this effect in 
models of RIS was, however, shown to be inconsistent 
with experiment^ for binding energies from 0.1 to 1.5 
eV and our results are consistent with these findings. 

The formation and binding energies of single Ni and 
Cr solutes to (100) self-interstitial dumbbells at Inn and 
2nn separation is presented in TABLE IXVII In both the 
fct afml and afmD states Ni is most stable in the ten- 
sile 2nn site, collinear with the dumbbell axis, exhibiting 
a binding up to approximately 0.1 eV. These are also 
the most stable configurations found here for a single 
Ni solute interacting with an overcoordinated defect, al- 
though in the fct afmD state the FeNi [001] dumbbell 
is close in energy. Although small, this binding allows 
Ni to act as a weak trap that inhibits interstitial diffu- 
sion and is therefore worth consideration in models of 
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Config. fee afmD fct afmD fct afml 

Ef Eb El E-b El _Eb 



[001 
[001 
[001 
[001 
[001 
[001 



[100 
[100 

[oio; 
[loo; 
[oio; 

[100 
[010 
[100 

[loo; 



[001 
[001 
[001 
[001 
[001 
[001 



[100 

[loo; 
[oio; 
[loo; 
[oio; 

[100 
[010 

[loo; 
[loo; 



-Ni/la 2.831 -0.074 3.286 -0.006 3.783 -0.001 
■Ni/lb 2.887 -0.129 3.249 0.030 N/A 

-Ni/lc 2.949 -0.191 3.357 -0.078 3.783 -0.001 
-Ni/2a 3.313 -0.034 3.861 -0.079 

■Ni/2b 3.291 -0.011 3.723 0.059 

-Ni/2c 3.209 0.070 N/A 



-Ni/la 
■Ni/lb 
■Ni/lb 
-Ni/lc 
-Ni/lc 
-Ni/2a 
-Ni/2a 
■Ni/2b 
-Ni/2c 



3.484 
3.373 
3.453 
3.379 
3.402 
3.279 
3.376 
3.404 
3.382 



-0.085 
0.027 
-0.053 
0.021 
-0.003 
0.121 
0.024 
-0.004 
0.017 



3.769 -0.071 

N/A 

N/A 
3.710 -0.012 
3.726 -0.028 
3.598 0.100 
3.747 -0.049 
3.692 0.006 

N/A 



Cr/la 2.953 -0.056 3.428 0.035 3.614 0.048 
Cr/lb 2.859 0.038 3.274 0.189 N/A 

■Cr/lc 2.843 0.053 3.270 0.193 3.547 0.115 

Cr/2a 3.500 -0.037 3.783 -0.121 

Cr/2b 3.394 0.069 3.670 -0.008 

-Cr/2c 3.467 -0.004 N/A 



Cr/la 
Cr/lb 
Cr/lb 
■Cr/lc 
■Cr/lc 
Cr/2a 
Cr/2a 
Cr/2b 
■Cr/2c 



3.407 
3.471 
3.492 
3.529 
3.461 
3.570 
3.552 
3.522 
3.598 



0.176 
0.112 
0.091 
0.054 
0.122 
0.013 
0.031 
0.061 
-0.015 



3.527 0.051 

N/A 

N/A 
3.577 0.001 
3.575 0.003 
3.526 0.052 
3.608 -0.030 
3.569 0.009 

N/A 



TABLE XVI: Formation energies, Ei, and, binding energies, 
-Bb, in eV for (100) defect to substitutional Ni and Cr solutes 
in Fe. Configuration labelling is as defined in FIG. 2(b) 



microstructure evolution and RIS, especially at higher 
Ni concentrations. Finally, it is worth noting that Ni is 
most strongly repelled from compressive sites at Inn to 
the self-interstitial dumbbell. 

For Cr there is less agreement between the two mag- 
netic states. In fct afmD the FeCr and CrFe [001] mixed 
dumbbells are the most stable, closely followed by Cr 
in compressive sites at Inn to a self-interstitial dumb- 
bell. The same configurations are found to be among the 
most stable in the fct afml state although others are at 
least as stable with no discernible preference for tensile 
or compressive sites. Despite these differences our dilute 
results show that, on average, Cr exhibits positive bind- 
ing to (100) dumbbells in mixed or neighbouring sites, 
where like Ni it would act as a weak trap for interstitial 
diffusion. 

Overall the interactions of Ni and Cr solutes with point 
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FIG. 8: Comparison of solute properties in austenitic and fer- 
ritic Fe. The dataset consists of the substitutional formation 
energy, solute- vacancy binding energies and solute- interstitial 
binding energies for the most stable self-interstitial, in eV. 
The black line is included to indicate an exact agreement. 



defects in austenitic Fe are consistent with them being 
moderately oversized and undersized solutes respectively. 
Ni binds to the vacancy and is generally repelled from 
the mixed and compressive sites of (100) dumbbells but 
shows positive binding in the tensile sites. Cr gener- 
ally exhibits the opposite tendencies. Such a conclusion 
should, however, be contrasted with the lack of consistent 
behaviour for single solutes in defect free Fe, as discussed 
earlier. Indeed, experimental results®^ find Ni to be un- 
dersized and Cr oversized in an austenitic Fe-Cr-Ni alloy, 
consistent with the pure elements. This may account for 
the earlier attempts to suggest a positive binding of Ni 
to interstitial defects. 

As a final means to draw conclusions from the data 
we have plotted our results for the solute substitution 
energy, solute-vacancy binding energies and binding en- 
ergies of solutes to the most stable interstitial dumbbells 
in the mixed positions in austenitic Fe against the cor- 
responding values for bcc Fe in FIG. [51 Ab initio data 
for bcc Fe was either as published previousl y ^"^' ^^ ' ^° or 
calculated here with similar settings. 

It is immediately clear that there is no discernible cor- 
relation between the two data sets, which is in com- 
plete contrast to the strong similarity found for pure Fe. 
The (vertical) spread of fee data also illustrates just how 
strong the effects of symmetry breaking within and be- 
tween the two afm datasets can be, especially when com- 
pared to the average binding energies. We suggest that 
this result simply emphasises just how sensitive the in- 
teractions of solutes and point defects are to their local 
magnetic environment in the Fe-Cr-Ni system in addition 
to the usual volume-elastic effects. The importance of de- 
veloping a deeper understanding of these effects cannot 
be overemphasised for the modelling of austenitic alloys. 
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VI. CONCLUSIONS 

We have carried out an extensive series of first princi- 
ples calculations to determine the energetics of austenitic 
steels. We first investigated a large set of possible refer- 
ence states for austenite at OK and found the fct afml 
and afniD states to be the most suitable, highlighting 
the problems with the unstable equilibrium fm-HS states 
in the process. It would be incorrect to associate any 
defect property of a paramagnetic material with a sin- 
gle microstate. However, by sampling various reference 
configurations we have been able to provide estimates 
of defect energies typically to within a few tenths of 1 
eV. The uncertainties associated with choice of refer- 
ence state in this work are in addition to the normal 
discrepancy between ab initio and experimental quanti- 
ties arising from, on one hand, choice of pseudopoten- 
tial, exchange-correlation potential, neglect of zero-point 
energy, and on the other from finite temperature and 
experimental errors. However, with considerably more 
computational effort, our dataset for austenitic materials 
has only moderately larger uncertainties than previous 
ones for ferritic steels. 

The main results and predictions of the paper can be 
summarised as follows: 

• The vacancy formation energy in austenite (at OK) 
is between 1.8 and 1.95 eV. Divacancy binding is 
rather weak at around 0.1 eV, suggesting that the 
nucleation of voids will face a nucleation barrier at 
elevated temperatures. There is, however, a rapid 
increase in total binding energy with cluster size 
e.g. 2.5 eV for a six vacancy octahedral cluster. 
Stacking fault tetrahedra exhibit similarly strong 
binding to the proto- voids (e.g. octahedral cluster) 
considered here and we therefore cannot distinguish 
between these two on energetic grounds. 

• The most stable self-interstitial is the (100) dumb- 
bell, consistent with other fee metals, with a for- 
mation energy of between 3.2 and 3.6 eV. These 
dumbbells aggregate strongly to form small defect 
clusters in (100) planes and although not consid- 
ered here, pair binding also suggests that clusters 
in (111) planes would also exhibit high stability. 

• Pair bond models for the total binding of defect 
pairs and clusters agree well with the data and are 
able to capture and highlight general trends. 

• Our results show that Ni and Cr do not strongly 
attract each other: binding energies of at most 0.1 
eV were found for Ni-Ni pairs. The binding en- 
ergies of solute clusters were consistent with pair 
interactions and fits to a pair bond model showed 
good agreement. Based on these results we sug- 
gest that short range ordering effects in austenitic 
alloys at operating temperatures typical of nuclear 
applications will be small generally but with some 
tendency to form Ni rich clusters. 



• We find that Ni binds to vacancies but by no more 
than 0.1 eV and can therefore act as a weak trap 
for vacancy migration, which may be of importance 
for void formation and RIS in concentrated alloys. 
The Cr-vacancy interaction is generally repulsive 
suggesting that its effect on microstructure evolu- 
tion and creep is negligible but concentration de- 
pendent effects cannot be ruled out in this study. 
Barrier heights for vacancy mediated migration are 
lowest for Cr, followed by Fe and finally Ni in the 
dilute case. This is consistent with the preferential 
diffusion necessary to explain the effects of RIS in 
austenitic alloys by the vacancy mechanism of the 
inverse Kirkendall effect. 

• Cr was found to bind to mixed interstitial defects 
whereas Ni is generally repelled. The preferential 
association of a particular solute with the inter- 
stitial flux under irradiation is considered an im- 
portant factor in RIS and our findings for Cr are 
worthy of consideration in the modelling of such 
effects. Substitutional Ni and Cr bind to (100) 
self-interstitial dumbbells at nearest neighbour sites 
and therefore act as weak traps for their migration. 

• Our results in austenitic Fe were compared with 
equivalent results in ferritic Fe. A strong similar- 
ity was found between these two datasets for point 
defect formation and binding energies in pure Fe. 
In contrast, no correlation was observed for Ni and 
Cr solute interactions with point defects. 

• We performed tests for the presence and infiuence 
of non-collinear magnetism in a subset of configura- 
tions. The tests either proved negative or resulted 
in marginal changes to the energy of the system. It 
should be borne in mind that such tests were useful 
to perform but are by no means exhaustive. 

In very broad terms, we have found that austenitic 
Fe behaves similarly to other fee metals, with (100) in- 
terstitials clustering to form proto-dislocation loops and 
vacancies clustering to form sessile stacking fault tetrahe- 
dra and voids. This normality is in contrast to the often 
anomalous behaviour of bcc Fe. The interactions of Ni 
and Cr with point defects are consistent with them being 
modestly oversized and undersized defects respectively 
despite the experimentally observed size factors and pure 
element data showing the opposite result. Migration of 
Cr through the lattice by both vacancy and interstitial 
mechanisms is enhanced relative to Fe self-diffusion. By 
contrast Ni migration is slow. Pure fee Fe is stable only at 
high temperature, but the similarity to other fee metals 
strongly suggests that FeCrNi steels stabilised at lower 
temperature by alloying will also behave normally. 

Of course, to make these predictions more quantitative 
it will be necessary to run further calculation using a 
method suited to larger systems and longer timescales, 
such as molecular dynamics or kinetic Monte Carlo. One 
of the main benefits of this study is to provide a large 
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database of configurations against which such models can 
be parametrised. 



tion. 
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